Measurement of Dissipation of a Three-Level rf SQUID Qubit 
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Abstract 

The dissipation-induced relaxation (Ti) time of a macroscopic quantum system - a A- type three-level rf SQUID flux qubit 
weakly coupled to control and readout circuitry (CRC) - is investigated via time-domain measurement. The measured interwell 
relaxation time of the qubit's first excited state, Ti — 3.45 ± 0.06 /is, corresponds to an effective damping resistance of the flux 
qubit R = 1.6 ± 0.1 MQ which is much lower than the intrinsic quasiparticle resistance of the Josephson tunnel junction. An 
analysis of the system shows that although the CRC is very weakly coupled to the qubit it is the primary source of damping. 
This type of damping can be significantly reduced by the use of more sophisticated circuit design to allow coherent manipulation 
of qubit states. 



PACS numbers: 85.25.Dq, 85.25.Cp, 03.65.Yz, 03.67.Pp 



The superconducting flux qubit, being immune to 
charge fluctuation, has attracted much attention in the 
area of solid-state quantum computation. The flux qubit, 
consisting of a superconducting loop interrupted by ei- 
ther one (rf SQUID) jlj] or three Josephson junctions 
(persistent current qubit or simply PC qubit) usually 
is a double well potential when a magnetic flux close to 
one half of a flux quantum is applied. The lowest states 
of each well serve as the two qubit states corresponding 
to macroscopic currents circulating around the loop with 
clockwise and anti-clockwise directions respectively and 
thus can be discriminated by a dc SQUID magnetometer. 
Quantum superposition of distinct macroscopic states, a 
prerequisite for quantum computation, was observed in 
both types of flux qubits at about the same time a couple 
of years ago Since then rapid development in PC 

qubit has led to the demonstration of coherent manipu- 
lations of a single qubit 0, entanglement between a PC 
qubit and its readout magnetometer and between two 
PC qubits (^]. Furthermore, the relaxation time (Ti) and 
dephasing time {T2) of the PC qubit have been measured 
directly using time-domain techniques providing invalu- 
able information about mechanisms of and solutions to 
decoherence 0,0,IEI3- contrast, despite intense ef- 
forts similar measurement have not been possible in the 
rf SQUID qubit due to shorter decoherence times in the 
particular systems measured. It is essential to study the 
sources of this decoherence to understand if it is funda- 
mental to this type of qubit or can be overcome by more 
sophisticated design. One very important element of this 
is the dissipation due to the coupling of the qubit to the 
environment. To the best of our knowledge the dissipa- 
tion rate of the rf SQUID qubit was inferred previously 
from escape probability distribution measurements which 
can provide only an order of magnitude estimate of the 
damping resistance i? 0, 0, Furthermore, the re- 

sults have led to more questions than answers because 
of the lack of an accurate knowledge about crucial sam- 



ple parameters and the use of an effective temperature 
that is more than ten times of the bath temperature to 
extract R from the data 0,^3- Therefore, a more di- 
rect and quantitative measurement of dissipation in rf 
SQUID qubits is needed in order to clarify the origin and 
the limit of the decoherence in the system. In this Letter, 
we report results of a relaxation time measurement of an 
rf SQUID qubit using time-resolved techniques. In our 
experiment all parameters that enter the qubit's Hamilto- 
nian are obtained from independent measurements which 
enables us to determine the damping resistance of the 
qubit with significantly improved accuracy. The result 
indicates that although the inductive coupling between 
the qubit and its control and readout circuits (CRC) is 
rather weak by conventional standards, it nevertheless is 
the dominant source of dissipation inducing relaxation 
from excited states. Hence, great care must be taken 
in the design of rf SQUID qubits and associated CRC to 
limit dissipation to levels allowing coherent manipulation 
of qubit states ■ 

The qubit used in our experiment is a variable bar- 
rier rf SQUID in which the single Josephson junction in 
an ordinary rf SQUID is replaced by a low inductance 
dc SQUID as shown in Fig. 1(a). The two-dimensional 
(2D) potential energy surface of such a variable barrier 
rf SQUID is 

U (0, (j)dc) = ($o/i) [{(1^ - 4>xf /2 + g{<t>dc - (t>xdcf/'2 

— /3o cos TT^dc cos 27r(/) + 6P sin iripdc sin 2tt(I>] (1) 

Here, $0 = h/2e is the flux quantum; g = L/21 is the 
ratio of the inductances of the rf SQUID and dc SQUID; 
/3o = 27rL/c/$o, 5f3 = 2ttL {1,2 - Id) /$o, 0^ (0^dc) and 
4' {4'dc) are the flux applied to and the net flux enclosed 
in the rf SQUID (dc SQUID) in units of $„• A plot of 
the qubit's first four energy levels is shown in Fig. 2(a). 
Notice that 0^ sets the energy bias e between the two 
wells while (f>xdc determines the tunnel splitting A. Thus 
the energy level structure of the qubit (e.g., e and A) can 
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FIG. 1: Scanning microscope image of the rf SQUID qubit. 
Inset: Equivalent circuit of the qubit inductively coupled to 
control (flux bias) and readout (magnetometer) circuits. For 
w/27r > 5 GHz « F„> ~ 
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FIG. 2: (color online) (a) Enrgy of the first four eigenstates 
of the qubit. Dashed line is the top of the potential barrier 
between the wells. Vertical arrows indicating the position 
of flux bias where the level spacing equal to 16.0 GHz. (b) 
Measured occupation probability of the \l) state as a function 
of flux bias without (empty circles) and with (solid triangles) 
16.0 GHz contonuous microwave irradiation. 

be varied in situ by adjusting (jj^ and (pxdc- The qubit 
is inductively coupled to a hysteretic dc SQUID magne- 
tometer and the state of the qubit can be determined by 
measuring the flux-dependent switching current of the 
magnetometer 0. 

In our experiment the rf SQUID is configured as a A- 
type three-level qubit with each potential well having 
only one level as the and 1 logic states and an aux- 
iliary level just above the potential barrier as shown in 
Fig. 2(a). Compared to a flux qubit in the usual two- 
level configuration the A-type three-level qubit has sig- 
nificantly faster gates and much lower intrinsic error rate 
[l^ . Moreover, since gate operations in such a three-level 
qubit are via the auxiliary level, it does not involve the 
tunnel splitting energy which depends exponentially on 
the potential barrier [l^- For these reasons the three- 
level flux qubit puts relatively less-strict constraints on 



sample fabrication. For e ^ A, the two lowest levels of 
the qubit are localized, corresponding to currents circu- 
lating in opposite directions in the qubit loop. For con- 
venience, we denote the state localized in the left (right) 
well by \l) {\r)). The energy difference between the |^) 
and |r) states can be varied continuously by sweeping 
the external flux . Microwave radiation with frequency 
Lo matching the energy difference E2 — Eq between the 
short-lived auxihary state |2) and ground state (e.g., |r) 
for e > 0) induces transitions between them. In general, 
the qubit can relax to the ground state from the auxiliary 
state via two pathways: |2) \l) \r) and |2) \r), as- 
suming s > 0. The probability of finding the qubit in the 
first excited state depends on various transition rates (in- 
cluding both stimulated absorption/emission and spon- 
taneous emission) between the three states involved. The 
lifetime of first excited state, Ti , is proportional to damp- 
ing resistance, R = 1/Re[y(a>)], where Y is the admit- 
tance seen by the qubit and u; is the transition frequency 
of the relevant radiative decay process 

[il,E3- There- 
fore, dissipation in an rf SQUID qubit can be investi- 
gated quantitatively via direct measurement of Ti be- 
tween qubit states. 

The sample was fabricated using a self-aligned Nb tri- 
layer process The size of the rf SQUID loop is 

50 X 100 /im^. In order to extract the qubit damping 
resistance from the measured Ti time one must have an 
accurate knowledge of key SQUID parameters, such as 
total junction capacitance C, loop inductance L, and 
/3o and Sp. In our experiment, the loop inductance 
L = 205 ± 5 pH and g — 17.0 ± 0.6 are estimated using a 
3D inductance calculation program, Pq = 3.70 ±0.02 and 
(5/3 < 0.05 were determined from the maximum size of 
the measured hysteretic 4> {4>x) loop (taking into account 
the effect of macroscopic quantum tunneling) and <p{<px) 
at (fixdc = 1/2, respectively [l^|20|. The total capaci- 
tance of the junctions in the rf SQUID, C = 65 ±2 fF, is 
obtained from microwave spectroscopy which agrees very 
well with the value determined from the total size of the 
junctions (1.3 /im^) and the specific capacitance of the 
Nb/AlOx/Nb trilayer (45 fF/^xm^). In addition, the crit- 
ical current density inferred from (3q, L, and the junction 
size, Jc — 460 A/cm^, agrees very well with that mea- 
sured directly from co-fabricated large junctions. The 
mutual inductances between the magnetometer and the 
qubit is Mm — 3.3 pH and that between qubit and flux 
bias line is = 1.0 pH, respectively. The sample is 
mounted in a oxygen-free copper cell thermally anchored 
to the mixing chamber of a dilution refrigerator. A super- 
conducting shield at 0.5 K, a cryoperm shield at 4.2 
K, and a ^-metal shield at room temperature are used 
to reduce flux bias fluctuations from the ambient mag- 
netic field. All leads to the sample cell are filtered with 
electromagnetic interference filters at room temperature, 
low-pass RC filters at 1.4 K, and microwave filters at mix- 
ing chamber temperature. A cryogenic coaxial microwave 
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FIG. 3: (color online). Excitation probability of qubit vs. 
microwave power at <j)x ~ 0.501 and uj = 2tt ■ 16.0 GHz. The 
solid line is the best fit to Eq. Q, where Pmw is measured 
at the top of cryostat. Inset: The potential (projected to the 
E-<j) plane), energy levels (dashed lines) and wave functions 
(solid lines) of the first three eigenstates. The arrows indicate 
some of the transition pathways. 



cable, with 30 dB and 20 dB attenuators thermally an- 
chored to the 1 K plate and mixing chamber respectively, 
couples microwaves to the sample. Battery-powered low- 
noise preamplifiers are used to monitor the bias current 
and voltage of the magnetometer. All ac-powered instru- 
ments are optically isolated from the battery-powered cir- 
cuits. Diagnostic tests using low critical current junctions 
(/c 1 /.tA) indicate that extrinsic noise is negligible at 
30 mK. 

We first measure the excitation probability, ni, of the 
state \l) vs. flux bias at constant microwave frequency 
and different A^dc to determine the energy level structure 
of the qubit [ij . By setting the amplitude and duration 
of bias current pulse properly the magnetometer either 
switches to finite voltage or stays at zero voltage depend- 
ing on the qubit state being |/) or |r). The result of the 
spectroscopy measurement is also used to select the val- 
ues of (j)x and (pxdc for time-resolved measurements. The 
result is shown in Fig. 2. For Pmw = 0, when s is 
swept from negative to positive there is a step-like tran- 
sition from n/ « 1 to n/ « indicating the change of 
ground state from \l) to \r). When irradiated by con- 
tinuous microwave (CW) a peak (dip) in n; appeared 
at ~ 1/2 + (— ) 0.0010, corresponding to microwave 
induced excitation from the ground state to the second 
excited state. The excitation are indicated by arrows in 
the energy level diagram of the qubit (Fig 2 inset), cal- 
culated using the independently determined sample pa- 
rameters and the full 2D potential When the value 
of (t>xdc is varied, while keeping uj constant, the position 
of microwave induced peaks shifted as expected from the 
calculated level diagram. 

Next, we measure n; as a function of CW power Pmw at 
(j)^^^ = 0.3868 where (Ea - EB.)/h = 27r ■ 16.0 GHz. Fig. 
3 shows the peak height as a function of the microwave 
power Pmw measured at the top of the cryostat. Since the 
microwave used in our experiment is cither CW or pulses 
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FIG. 4: (color online). Excitation probability of the first 
excited state vs. pulse width at three different levels of mi- 
crowave power. The solid lines are best fits to Eq. Q. The 
resulting T\ is independent of the microwave power applied. 
Inset: Signals employed in time-domain measurement of qubit 
excitation probability. Top: Amplitude of microwave applied 
to the qubit. Middle: Bias current of the magnetometer. The 
switching current of the ground (first excited state) state is 
greater (less) than the maximum bias current. Bottom: Volt- 
age across the magnetometer. 

with duration much longer than the decoherence time, 
the dynamics of the qubit is incoherent. Hence, the time 
evolution of the qubit state under microwave irradiation 
can be described by the following master equation 

dno/dt = -Tno + 7ioni -f (F + 720) ^2, (2a) 

dni/dt = -7io?^i + 72lf^2, (2b) 

9712/c'i = Fno - (r + 721 +72o)".2, (2c) 

where, no (i), ni (t), and n2 (t) are the normalized pop- 
ulations of the ground state, the first excited state, and 
the second excited (the auxiliary) state, F is the rate 
of stimulated transitions between the ground state and 
the auxiliary state, and 7^ is the spontaneous decay rate 
from \i) to \i < i) — 0,1,2). Note, in our experi- 
ment we have |; = 0) (|r = 0)) for 0^ < 1/2 {(j)^ > 1/2) 
while |a = 2) regardless the value of <j)x- For weak mi- 
crowave fields one has F = aPmw according to Fermi's 
golden rule, where a is a coupling constant dependent of 
the microwave circuit used. Under irradiation of CW, the 
three-level qubit is in a steady state, dni=i^2.3/dt — 0, 
and Eq. (2) can be solved exactly to give 

p 

^ mw 

(1 -I- 2710/721) Pm«; + (710/721) (720 +721) /a' 

(3) 

Fitting the measured rii vs. Pmw to Eq. Q yields 
710/721 = 0.072 ± 0.040 and (720 + 721) /a = 4.0 ± 0.8 
mW. This value of 710 /721 corresponds to a ratio between 
the two relevant matrix elements |0io/</'2i| = 0.26 ±0.05 
which agrees very well with \(j)io/4>2i \ = 0.254 calculated 
using the independently determined qubit parameters. 

Although the steady state measurement provided use- 
ful information about the system one cannot extract 
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Ti = from it alone. To obtain Ti we used a time- 
domain technique. We start each measurement cycle by 
applying a microwave pulse of duration tp. A bias cur- 
rent pulse is then applied to the magnetometer to read- 
out the state of qubit. The flux bias is kept constant 
at — 0.5010. The time evolution of rt; is obtained 
by varying tp and repeating the cycle more than 1000 
times for each value oi tp. A few milliseconds of delay 
was inserted between successive cycles to allow the qubit 
to relax to the ground state before each new measure- 
ment cycle. Fig. 4 shows ni{tp) for three different mi- 
crowave power levels. Since F/ (720 + 721) = Pmw/{4:-0 
mW) < 0.1, we have F <C 720 + 721 ■ Under this con- 
dition, 712 (t) rapidly saturates to a constant n2oo, which 
depends on P^w, at t > to = 1/ {2r + 721 +720) < 7iV 
and (|2bp can be solved analytically, by replacing 712 (t) 
with n2oo, to give: 

ni (t) = n2oo(7io/72i)"' (l - e"^"*) . (4) 

We verified the validity of the approximate solution 
by solving Eq. (2) numerically. The result is in excellent 
agreement with at t > to. Since 710/721 is known 
from the steady state measurement, fitting the measured 
ni{t) to Q gives Ti = -f^J = 3.45 ±0.06 ^s. In addition, 
the results of steady state and time dependent excitation 
measurements together yield to ~ 70 ns and a « 3.0 
mW^^/iS~^ confirming the applicability of iQ. 

It is well known that the damping resistance's contri- 
bution to Ti time of an rf SQUID qubit is proportional 
to i? Hi 

Tr' = 710 = i27T/h){Ei - Eo){Rq/R) 

{l + coth[(£;i -So)/2fcBT]}, (5) 

where Rq = h/Ae^ is the resistance quantum. From 
(HI, using the value of <^io = 1.0 x 10~^ calculated 
from the independently determined qubit parameters, we 
have R = 1.6 ± 0.1 Mfi, which is more than 10^ times 
lower than the measured quasiparticle resistance of co- 
fabricated junctions Hence, we conclude that the 
contribution of quasiparticles to qubit dissipation is neg- 
ligible. On the other hand, weak but finite inductive cou- 
pling to the control and readout circuits result in a damp- 
ing resistance R ~ {Re[Fm {t^w)+Yx (t^io)]}~^j where Ym 
and Yx are the admittances of the magnetometer and the 
flxix bias circuits seen by the qubit and loio = {Ei~Eo)/ h 
is the transition frequency between the |^) and \r) states, 
respectively. For the CRC used in our experiment which 
is illustrated in Fig. 1(b) 

~ Re(y,„ + r,) 

~ Ro^[{M^/Lf{ALj/2E„,f + {Mx/Lf] (6) 

around ujiq, where ALj is the difference of Josephson 
inductances of the two junctions in the magnetometer 



and Lm is the sum of the geometric inductance 40 
pH) and the Josephson inductances of the magnetometer 
loop. Since the high frequency dampings (at to ^ ujio) of 
the flux bias and magnetometer's bias/measurement cir- 
cuits are essentially the same, we model them by a shunt- 
ing impedance i?o- Substituting the values of M^, M^, 
L, and Josephson inductances of the two junction evalu- 
ated at working point of the magnetometer (lb = 1.0 /lA, 
(j)m — 0.45) into (O yields Ro = 69 ± 5 ft, which agrees 
very well with i?o ~ 70 f2 derived from the measurement 
of escape rate of the magnetometer in phase-diffusion 
regime . The value of Ro obtained is typical of the 
high frequency impedance of transmission lines indicat- 
ing that the interwell relaxation of the qubit is mainly 
induced by its coupling to electromagnetic environment 
through the flux bias and readout circuits. In principle, 
changing the working point of the magnetometer and fur- 
ther reducing coupling between the qubit and CRC could 
result in longer coherence time. Unfortunately, it could 
not be done in our experiment since the magnetometer's 
maximum critical current, Ico = 9.50 fiA, is more than 
three times higher than the target and consequently the 
flux applied to the magnetometer {(pm) is bounded in a 
very tight range slightly below $0/2 to achieve nearly 
single-shot readout. On the other hand, according to 
Eq. jnj, contribution to the damping resistance from the 
integrated on-chip flux bias line is 1/Re(l^) ~ 2.9 MSI. 
The use of off-chip flux bias, as those employed in the 
experiments observing coherent oscillations in the persis- 
tent current qubits |j, H, ' would significantly increase 
1/Rc(y3;) but this approach would be very difficult to 
apply to a circuit containing many qubits required for 
practical quantum information processing. Currently, 
we are developing advanced designs for the qubit bias 
and readout circuits that are predicted to decrease their 
contributions to the qubit damping by several orders of 
magnitude. In addition, due to limitations in our mi- 
crowave coupling circuit, this rf SQUID qubit was con- 
figured as a magnetometer which makes it susceptible 
to ambient field fluctuations. This source of decoherence 
can be greatly suppressed by using qubits with gradiome- 
ter geometries. 

In summary, the dissipation of a Nb rf SQUID qubit, 
which is coupled inductively to the flux bias and readout 
circuits, is determined by measuring the interwell relax- 
ation time between the ground state and flrst excited 
state with the A-type three-level configuration. Time- 
domain measurement of the excitation probability of the 
first excited state yields Ti = 3.45 ± 0.06 /iS correspond- 
ing to a damping resistance R — 1.6 ± 0.1 Mfi for the 
qubit Analysis of the system indicates that the dominant 
sources of qubit dissipation are the fiux bias and magne- 
tometer readout circuits. Since this kind of dissipation- 
induced qubit decoherence can be greatly suppressed 
with more sophisticated designs we believe it does not 
impose a fundamental limit to this type of qubit. 
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